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Abstract Supermassive binary black holes (SBBHs) are 
a natural outcome of galaxy mergers. Here we show that 
low-frequency (/ < 10~ 6 Hz) quasi-periodic variability 
observed from cosmic blazar sources can provide sub- 
stantial inductive support for the presence of close (d ~ 
0.1 pc) SBBHs at their centers. It is argued on physical 
grounds that such close binary systems are likely to give 
rise to different (although not independent) periodicities 
in the radio, optical and X-ray/TeV regime, and, hence 
that detection of appropriate period ratios significantly 
corroborates the SBBH interpretation. This is illustrated 
for a binary model where optical longterm periodicity is 
related to accretion disk interactions, radio periodicity to 
Newtonian jet precession, and periodicities in the high 
energy bands to the orbital motion of the jet. We use 
the observed periodicities to constrain the properties for 
a sample of SBBH candidates including OJ 287 and AO 
0235+16, and discuss the results within the context of 
jet activity and binary evolution. 

Keywords Supermassive binary black holes, periodic- 
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1 Introduction 

Violent galaxy mergers are known to play a vital role 
in the cosmological evolution of galaxies and the growth 
of supermassive black holes (BHs). Since almost every 
bright galaxy seems to contain a black hole (BH), fre- 
quent formation of supermassive binary black holes is 
naturally expected to occur during cosmic time, e.g., 
see [1,2,3 . In cosmic bottom- up scenarios, for exam- 
ple, elliptical galaxies, such as the host galaxies of ra- 
dioloud AGNs, are usually thought to form when two 
spiral galaxies collide and merge. Interacting galaxies 
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are indeed observationally well-known and the obser- 
vation that the relative number of spiral to elliptical 
galaxies tends to increase in distant clusters [4], gives 
additional credit to the underlying evolutionary picture. 
Moreover, direct observational evidence for formation of 
wide (separation d> 1 pc) supermassive binary systems 
has been established recently based on Chandra obser- 
vations of the ultraluminous infrared galaxy NGC 6240 
[5] and VLBA observations of the radio galaxy 0402+379 
G . While the existence of wide SBBHs may thus be re- 
garded as observationally well grounded, the anticipated 
existence of close (d ~ 0.01 — 1 pc) SBBHs appears much 
more ambiguous. Indeed, there is little direct observa- 
tional evidence so far for a close, secondary supermas- 
sive BH in the nuclear region of Active Galaxies. On 
the other hand, dynamical friction and slingshot inter- 
actions with stars normally ensure that a binary system 
gets quickly closer, reaching separations d c ~ 1 pc in less 
than 10 s yr [lj. As the number of field stars with ideal 
impact parameters decreases with d, the binary evolu- 
tion will then stall at around d c (i.e., well above the 
separation at which gravitational radiation becomes im- 
portant), unless further angular momentum can be effi- 
ciently removed by other processes such as, for example, 
influx of gas or accretion disk interactions. SBBH sys- 
tems are thus likely to spend most of their time at sepa- 
rations d ~ (0.05 — 1) pc for BH masses of order 10 8 Mq. 
While the origin of X-shaped radio morphologies in some 
radio galaxies may provide some phenomenological evi- 
dence for a possible spin flip during the final merger stage 
[7], thus suggesting that at least some SBBHs may co- 
alesce within a Hubble time, it is theoretically still not 
yet fully understood today whether that can be achieved 
by a substantial fraction of SBBHs, e.g., see [i~1l8"ll9lll"0l 
ITT | [T2 l [T3] . As shown below, our analysis presented here 
provides further phenomenological support for relatively 
short merging timescales in AGNs. 
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2 Periodic variability in close SBBHs 

The presence of close SBBH systems has been repeat- 
edly invoked as plausible source for a number of ob- 
servational findings in blazar-type AGNs, ranging from 
misalignment and precession of jets to helical trajecto- 
ries and quasi-periodic variability, see [11,15 for recent 
reviews. As illustration of the latter, consider a simple 
SBBH model with a precessing jet: 

Interactions of the companion with the accretion disk 
around the primary BH may then naturally account for 
longterm optical periodicity with periods of the order of 
~ several years as observed in a number of blazars, 
e.g., see [16] for OJ 287 and also Table ([T]), at least in 
those cases where the disk provides a non- vanishing con- 
tribution to the observed optical spectral flux. Accord- 
ingly, we may derive an upper limit for the intrinsic Ke- 
plerian orbital period of the binary 

p, < ? p°p* (i) 

by assuming that the optical longterm periodicity is caused 
by the secondary BH crossing the disk around the pri- 
mary twice per orbital period. Note, however, that due 
to internal disc warping and/or disc precession, some de- 
viations from strict periodicity are likely to occur. 
The detection of helical jet paths in an increasing num- 
ber of blazar sources [T7lfT8l[T9l[20] suggests that quasi- 
periodic variability - especially in those energy bands 
dominated by the jet, e.g., radio, X- and 7-ray - may also 
naturally arise as a consequence of differential Doppler 
boosting S v (t) = 5{t) n S' v for a periodically changing 
viewing angle, where S(t) is the time-dependent Doppler 
factor, n > 3, and S v ' the spectral flux in the comoving 
frame [2"Tll2"2"] . For non-ballistic helical motion, classical 
travel time effects will then lead to a shortening of ob- 
servable periods P bs with respect to the real physical 
driving period P such that 
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where jb — (5 — 15) is the typical bulk flow Lorentz 
factor [2D] . Orbital motion and (Newtonian) disk preces- 
sion caused by tidally induced perturbations in the disk 
[2"3"I|2"4] belong to the most obvious driving sources for 
helical jet paths. If, as it is usually believed, the high 
energy emission is produced on small jet scales, it may 
be primarily modulated by the orbital SBBH motion, 
suggesting observable periods 
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as indeed observed [cf. Table dTJ)]. On the other hand, 
simple cooling arguments suggest that a significant part 
of the radio emission may originate on larger scales, where 
Newtonian jet precession can no longer be neglected. In 
general, the driving period for Newtonian jet precession 



is (at least) an order of magnitude higher than the orbital 
period Pk, i.e., P = P p > 10 Pk 20 , so that observable 
radio periods are expected to satisfy 

Pj b a s dio ^20P o r/7 6 2 - (4) 
Note, that if P p is rather small (say P p ~ 10 Pk), mod- 
erate bulk Lorentz factors are still sufficient to account 
for Pobs'° < P°bg as sometimes observed, cf. Table ([!]). 
Let us note, that a close SBBH interpretation for the 
origin of (some) quasi-periodic variability is certainly 
not the only possible explanation, as other (not mutu- 
ally exclusive!) origins (e.g., disk instabilities, orbiting 
disk hot spots) are conceivable as well. What makes the 
SBBH interpretation unique, however, is that it seems 
to allow corroboration of facts otherwise not possible: 
(1.) It is based on quite general arguments for (bottom- 
up) structure formation (galaxy mergers), (2.) it nat- 
urally accounts for helical jet trajectories observed in 
many sources, (3.) it can offer a reasonable solution to 
the problem of divergent central mass estimates derived 
from high energy emission models and host galaxy obser- 
vations [25], and (4.) it provides a coherent explanation 
for longterm periodic variability. In particular, as shown 
above, strong support can be provided by the detection 
of quasi-periodic variability on different timescales in dif- 
ferent energy bands. 



3 Application to individual sources 

3.1 AO 0235+16: 

Long-term monitoring (1975-2000) of this well-known 
and highly variable BL Lac object has shown evidence 
for a (5.7 ± 0.5) yr radio and a possible (2.95 ± 0.15) 
or (5.7 ± 0.5) yr optical periodicity 26,27 . Both find- 
ings have been interpreted within a close SBBH frame- 
work [28]l29j. cf. [30] for a discussion. However, during 
the latest radio to optical monitoring campaign in 2003- 
2005, no evidence for a major radio or optical outburst 
- extrapolated from previous observational results to oc- 
cur within the campaign - was found [3"l"]l3"2] . While the 
non-detection of a major radio outburst might be rel- 
atively easy accommodated in a precession-driven heli- 
cal jet model [28] by taking, for example, a moderate 
change of the jet bulk Lorentz or Doppler factor (de- 
crease) and/or the inner disk properties during an ac- 
tive stage fully into account, the non-detection of the 
expected optical outburst (if indeed associated with ac- 
cretion disk interactions, e.g., [28]), is somewhat more 
challenging for a consistent SBBH interpretation. A pos- 
sible way out - apart from the possibility of inaccurate 
previous periodicity results, i.e., the real optical period 
may actually be ~ 8 yr [32], which needs to be checked 
in detail - is to assume that, similar to the accretion- 
ejection connection in microquasars, an active (low-hard- 
type) source stage may be associated with enhanced in- 
trinsic (!) jet activity and a decrease in optical disk flux, 
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Table 1 Properties of a sample of blazar SBBH candidates, cf. Rieger et al., in preparation for more details, with masses in 
10 s Mq , and upper limits for the associated binary lifetime T gra v due to gravitational radiation in units of 10 s yr for typical 
mass ratios > 0.01. References given are for observed periods P b s identified. 



name 


redshift z 


periods P obs Ref. 


(m + M) 


Pk [yr] 


d/10 lb cm 


Tgrav 


Mkn 501 


0.034 


23.6 d (X-ray) 


43 




(2-7) 


(6-14) 


(2.5-6) 


5.50 






~ 23 d (TeV) 


44 












BL Lac 


0.069 


13.97 yr (optical) 


38 




(2-4) 


(13-26.1) 


(4.8-9.7) 


28.9 






~ 4 yr (radio) 


45 












ON 231 


0.102 


~ 13.6 yr (optical) 


46 




> 1 


(12.3-24.7) 


> 3.7 


79.4 






~ 3.8 yr (optical) 


47 










3C 273 


0.158 


13.65 yr (optical) 


48 




(6-10) 


(11.8-23.6) 


(6.5-12.3) 


3.55 






8.55 yr (radio) 


49 












OJ 287 


0.306 


11.86 yr (optical) 


37 




6.2 


(9.1-18.2) 


(5.5-8.8) 


1.68 






~ 12 yr (infrared) 


50 
















~ 1.66 yr (radio) 


42 
















~ 40 d (optical) 


39 












3C 66A 


0.444 


4.52 yr (optical) 


27 




> 1 


(3.1-6.3) 


> 1.5 


2.08 






65 d (optical) 


51 










AO 0235 


0.940 


2.95 yr (optical)? 


27 




> 1 


(1.5-3.1) 


> 0.95 


0.31 






8.2 (optical)? 


32 






(4.2-8.4) 


> 1.81 


4.48 






5.7 yr (radio) 


26 












3C 446 


1.404 


4.7 yr (optical) 


52 




6.5 


(1.9-3.9) 


(2.0-3.2) 


0.03 






5.8 yr (radio) 


53 













so that the corresponding disk contribution might be 
swamped by emission from the jet. If this is the case, 
then one might expect some orbital-driven optical vari- 
ability [cf. eq. (J3])] to occur on timescales ~ 20 d assuming 
7b, optical — 8 for the optical regime [33J . Clearly, exten- 
sive multiwavelength monitoring of this source around 
the next SBBH anticipated outburst and advanced anal- 
ysis methods will be important to clarify whether the 
latest anomaly may turn into a falsification of a sim- 
ple SBBH scenario, where optical longterm periodicity 
is related to disk crossing. We note that circumstantial 
evidence for very high Doppler factors ~ 100 (76 > 50) 
in AO 0235 have been reported in the literature, e.g., 
see [34] . which - if indeed true - would even allow for 
some orbital-driven intraday variability. A very interest- 
ing alternative explanation for the origin of some radio 
QPOs in AO 0235+16 (including its ~ 5.7 yr period- 
icity) has been proposed recently [35], suggesting that 
they might be related to periodic plasma injection into 
the jet driven by the p-mode oscillation of a thick inner 
disk (probably excited by a close SBBH) with a rela- 
tively high transition radius R tI ~ 10 3 r g (cf. however 
also the general arguments for much smaller transition 
radii in BL Lacs, [36] ) and associated intrinsic fundamen- 
tal frequency / ~ 0.18 yr" 1 ^/!^ • 10 3 r g )~ 3 / 2 (4.7 • 
1O 8 M /M B h)- As the QPOs are thought to arise due 
to periodic plasma injection, the model seems to im- 
ply that similar QPOs should be observable in different 
jet-dominated energy bands (i.e, not only in the radio), 
which may allow a straightforward test. To complicate 
matters, note however that due to the above noted travel- 
time effects the observed radio periods may not neces- 
sarily correspond to the intrinsic driving periods. 



3.2 OJ 287: 

Optical and infrared monitoring of this famous BL Lac 
object have shown strong evidence for a ~ 12 yr [371 
[38] longterm, and a possible ~ 40 d midterm periodic- 
ity [39] . with the 12 yr periodicity commonly interpreted 
as due to a close SBBH system, e.g., see [371IT61I40] . As- 
suming the 12 yr QPO to be caused by disk crossing 
implies Pk ~ 18.2 yr. Interestingly, orbital-driven helical 
jet motion [Eq. ([3])] then suggests observable midterm 
periodicity P ODS ~ 50 d for 7^ ~ 13 as derived from 
SED multiwavelength modelling of OJ 287 [41], which 
is well consistent with the observed timescale. If the ra- 
dio emission emerges from larger Newtonian-precession- 
modulated scales with P p ~ lOPfc, radio QPOs with 
-fobs ~ 1.5 yr might be expected assuming a similar bulk 
Lorentz factor for the radio regime. Again, there is sig- 
nificant evidence for such a periodicity in the radio data 
[42], suggesting that a close SBBH may indeed offer a 
powerful explanatory framework. 



4 Jet activity and evolution of SBBHs 

Table {1} presents properties derived for a sample of 
blazar SBBH candidates, where the observed periodic- 
ities have been used to estimate the last three columns: 
In all cases the gravitational lifetime of the binary is (a) 
significantly smaller than the Hubble time for the most 
likely range of mass ratios, a result that gives strong phe- 
nomenological support to the notion that SBBHs can in- 
deed coalesce, and (b) still large enough to satisfy the 
minimum source lifetime required for jet fuelling (e.g., 
t ~ 10 6 yr for the quasar 3C 273). The derived separa- 
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tions d are of order of the maximum size ~ 1000 r g 
(as given by the Toomre stability condition Q ~ 1) 
for a standard Shakura & Sunyaev disk around the pri- 
maryQ Accordingly, our results suggest that the candi- 
date sources may undergo early stages of binary-disk in- 
teraction, where the orbit of the secondary BH is still in- 
clined with respect to the circumprimary disc [55]. This 
may qualify our assumption that the secondary BH in- 
deed hits the disk around the primary twice per orbital 
period, so that the larger (upper bound) values for 
and d in Table ([I]) appear to be more realistic. Accretion 
disk interactions will lead to an accelerated evolution of 
the binary and can also excite the binary to eccentrici- 
ties e ~ 0.1 [55|fl2] . Gas accretion and binary-disc inter- 
actions are thus very likely to dominate the most critical 
binary evolution stage between d ~ 0.01 — 1 pc, support- 
ing the notion that SBBH systems are important but 
temporary feature of galaxy evolution, e.g., [9]. Binary- 
disc interactions can also provide a natural trigger for 
enhanced jet activity (e.g., time-dependent increase of 
accretion rate, cf. [37]) and may even lead to recurrent 
ejection of superluminal jet components. Moreover, dur- 
ing collision the disk gas will be perturbed, shocked, 
heated and ejected at the point of collision, which may 
lead to the possible formation of hot, optically thick out- 
flows from the disk with maximum luminosity up to the 
Eddington limit of the secondary [56] . 



5 Conclusion 

While each piece of observational evidence (e.g., helical 
jet trajectories, periodic variability etc.) may in principle 
allow a variety of alternative interpretations, we believe 
(cf. §2) that a strong cumulative (phenomenology-based) 
case can be built for the presence of close supermassive 
binary black holes at the centers of (at least some) ra- 
dioloud Active Galaxies. If so, then simple AGN uni- 
fication schemes may be fundamentally incomplete as, 
for example, activity cycles, jet structure and emission 
properties may be strongly affected by the presence of a 
secondary BH. Longterm monitoring of blazar sources, 
the use of astronomical plate archives and sophisticated 
analysis methods may thus be crucial for further in-depth 
assessment of the robustness of the observational basis. 
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